Instabilities of the Abrikosov-Suhl resonance 
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ABSTRACT: We consider the possibility that instabilities of the Abrikosov-Suhl reso- 
nance lead to new fixed point behavior of the Kondo effect in a lattice environment. In one 
scenario, a pairing component to the resonant scattering develops in the Kondo lattice, 
leading to an odd frequency superconductor. We discuss experiments that can discriminate 
between this picture and d-wave pairing, and its relationship to the non Fermi liquid fixed 
point of the overscreened Kondo model. 



Paper presented at the SCES conference, San Diego, 1993 in celebration of Harry Suhl's 
70th Birthday. 



Harry Suhl's birthday provides us with a welcome opportunity to pause and look back 
at early work on the Kondo effect in which he played a vital role. The Kondo or "s- 
d" model for magnetic ions in a metallic host dates back to work by Zener 1 ' 2 . In 1956, 
Kasuya 3 first cast the model into its modern second quantized form 

H = H c + J^Bj ■ Sj (<Jj = ^jaVaP^jp) (1) 

j 

Here H c = ^ ka^ka describes the conduction band and the exchange interaction is 
written in a tight binding form. In the early sixties Kondo, 4 building on Anderson's su- 
perexchange concept 5 ' 6 , studied the properties of the model with antiferromagnetic rather 
than the ferromagnetic interactions envisaged by Zener. The famous logarithmic correction 
to the electron scattering rate that he derived 

- = n t -K J 2 p \l + 2(Jp) Zn(-^)j (2) 
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where p is the conduction electron density of states and the impurity concentration, 
beautifully explained the resistance minimum of dilute magnetic alloys and definitively 
confirmed the indirect antiferromagnetic exchange between magnetic ions and conduction 
electrons. 

Suhl was one of the first to appreciate the many body implications of this logarithmic 
term. In the first of two papers published in 1965 7 he remarked "A divergence of this 
sort calls into question the stability of the Fermi surface" . Nowadays, we view the Kondo 
model in the language of scaling theory, and the logarithmic terms that prompted Suhl's 
remark are taken as an indication that the original conduction sea with unquenched local 
moments is an unstable fixed point. Perturbation theory tells us that spin fluctuations are 
"anti-screening" causing a flow of the superexchange to strong coupling as a function of 
energy, according to the perturbative beta function 

!^=/3(Jp), 0(x) = -2x 2 + 2x 3 + . . . (x«l) (3) 
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SuhPs early work on the impurity model showed that the logarithmic growth of the electron 
scattering amplitude leads ultimately to the development of an elastic resonant scattering 



understood to be a renormalized Friedel- Anderson resonance. 

Suhl's remark acquires a renewed significance in the light of the discovery of heavy 
fermion compounds. A priori, the presence of the Kondo logarithms in weak coupling tells 
us nothing about how more general Kondo models flow to strong coupling: this depends 
on the topology of the scaling flows. Experimentally, most heavy fermion compounds do 
indeed show features that demonstrate, that to some degree or another, their behavior 
is dominated by a flow to a Fermi liquid fixed point. However, almost all of them show 
magnetic or superconducting instabilities at temperatures that are substantial fractions of 
the Kondo temperature: despite the itinerant aspects to these phase transitions, they are 
in essence spin ordering processes, and require descriptions where the correlation of the 
local moments with the conduction electrons is an integral part of the ordering process. 

These considerations motivate us to consider possible instabilities of the Abrikosov 
Suhl (AS) resonance in more general Kondo models where the appearance of new relevant 
variables in the Hamiltonian divert the screening process into new basins of attraction. A 
simple example of this phenomenon is the impurity model with an additional screening 
channel coupled to the local moment, 



Here ct labels the screening channel. A separatrix at A = 1 divides the Hamiltonian flows 
into two distinct Fermi liquid basins of attraction where the AS resonance resides in one 
channel, but is absent from the other. At A = 1 the AS resonance becomes unstable, 
and the model scales to a non Fermi liquid (NFL) quantum critical point, developing a 
unique localized real fermion mode with a fractional entropy. 12 (Fig. la.). There has been 



center at the Fermi energy. The "Abrikosov- Suhl" resonance that he predicted 7 



is now 




(4) 
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much recent interest in this fixed point in connection with the possibility of a quadrupolar 
Kondo effect in heavy fermion systems. 1 ^ 

Does the lattice play a similarly relevant role in the Kondo scaling process? In the 
lattice, crossing a separatrix between two basins of attraction would imply a real phase 
transition(Fig. lb.) associated directly with the Kondo effect. One interesting possibility 
is an instability in the Abrikosov-Suhl resonance that leads it to develop an anomalous scat- 
tering component in the triplet channel, producing an odd frequency pairing component 
in the electron self-energies. 



This type of "odd frequency" triplet pairing was envisaged by Berezinskii 14 ' . Calculations 
that we now outline 16 indicate that such a state is stable, sharing features in common with 
the fixed point of the overscreened Kondo model, notably the development of neutral 
fermionic modes at the Fermi energy that decouple from the spin and charge degrees of 
freedom. The state that is formed retains certain features reminiscent of a superconductor 
containing a line of gap zeros, making it a conceivable alternative to the d-wave theory of 
heavy fermion superconductivity. 

A key step in our work is to bypass the difficulties of gauge theory approaches to the 
Kondo lattice model, replacing the commonly used Abrikosov pseudo-fermion representa- 
tion of spins with a more ancient, Majorana representation that avoids the constraints. 
This method, in a different guise, was used by Spencer and Doniach 17 in their early work 
on the Kondo model under the name of "Drone fermions" , but has since been neglected. 
The essence of the method is to generalize the fermionic properties of Pauli matrices to 
a lattice. Recall that Pauli matrices anticommute {a a , 05} = 25 a b which implies they are 
real or "Majorana" (<jt = a) fermions. Their Fermi statistics guarantee that the "spin" 
operator S = —^axa is a faithful representation of a spin-1/2. We generalize this property 





(5) 
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to the lattice by introducing a three component Majorana fermion fji at each site i, which 
satisfies the algebra {t]f , rjj} = Sijd ab , (a, b = 1,2,3). The corresponding spin operators 
are then 

$j = ~\m x m ( 6 ) 
In terms of Majorana fermions, the Kondo exchange interaction can be written as 

H int [j] = -jV^Vj (7) 

where 

Vj = (j}A = -^-m (8) 

is a charge e spinor formed between the conduction electron and local moment. This 
form suggests the possibility that electrons and local moments will condense to develop a 
vacuum expectation value 

VjA-fm^' i (9) 



V 3lJ V<0|%| 

Defects of a charge e spinor correspond to tt changes in the order parameter phase and 
carry the same flux quantum as a charge 2e scalar. Spinor order of this type generates a 
non-trivial correlation between the spin of the local moments and the conduction electron 
pair degrees of freedom 

(T a (x)SP(x)) = gM£(x) (a, (3 = 1, 2, 3) (10) 

Here r(x) is the conduction electron "isospin" : z components describe the number density, 
r 3 = \{p{ x ) ~ 1)? an d transverse components describe the pairing t+(x) = tp^ (x)tp^(x). 
The quantity g ~ Vj / J 2 defines the magnitude of the order parameter and M_ is an 
orthogonal matrix whose orientation is set by the components of V. 

On a bipartite lattice, the lowest energy stable mean field solution is produced by a 
staggered pairing field of the form 

v e i{Q-Rili)V z (11) 
V2 



where Q = (tt, it, it) and Z is a unit spinor. The mean field Hamiltonian describes a mixing 
between local moments and conduction electrons, where the hybridization is a spinor: 

H mft = E Vfc + ^ [^V* ' + H ' C -] } + (12) 

(Here the staggered phase has been gauge transformed to the conduction electrons by the 
replacement — > = ^^qj^}- When a conduction electron hybridizes with the zero 
energy Majorana modes, it develops a self-energy component proportional to 1/u. Since 
the Majorana fermion is neutral, the scattered fermion can emerge as either electron, or 
hole (see below), thereby developing an odd frequency pairing term in the conduction 
electron self energy. 

Diagrammatically, this process is represented as 

V2 i 

> — • • > = — 1 

y2 ' ( 13 ) 

— >— <: = 2^ 

where a dotted line indicates the intermediate resonance and we have used a Nambu 

notation to denote anomalous pairing components. For example, if Z = ( . J then the 



down electrons experience this resonant pairing, with self energy 



V 2 

E^) = —V (14) 

where V = \ [l + 7jJ is a projection operator that projects out a Majorana component of 
the down conduction sea: the remaining half does not couple to the odd-frequency triplet 
pairing field and forms a novel band of decoupled gapless excitations. If we decompose the 
conduction electron spinor into constituent "Majorana" components 

1 



1>l + i1> v az (15) 
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the zeroth component, ip*Z = — "^[V^ ~ V^^] is decoupled from the local moments. The 
other three components of the conduction bands develop a gap A g ~ V 2 /D, (where D 
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is the bandwidth) that decouples them from the gapless mode, rendering it both neutral 
and spinless. The gapless excitations consequently have vanishing spin/charge coherence 
factors at the Fermi energy which vanish linearly with energy 



where /i is the chemical potential and uo^ the quasiparticle energy. Linear coherence factors 
lead to power laws in the nuclear magnetic relaxation 

Ti 

that coexist with a linear specific heat capacity. Unlike a d-wave superconductor with 
lines of zeroes, this T 3 NMR relaxation rate will persist even when the linear specific 
heat is large. In a dirty d-wave superconductor we expect a Korringa relaxation in the 
superconducting state once it develops a linear specific heat. Fig. 2. contrasts the linear 
specific heat and NMR relaxation rate for the simple mean field theory outlined here, 
demonstrating these features. 

As in all mean field theories, issues of stability and fluctuations are of paramount 
importance. Uniform odd- frequency states appear in general to be unstable. 16 Here the 
staggered phase stabilizes the state, producing a finite Meissner stiffness. Since there are 
no awkward gauge modes, fluctuations about the odd-w state are similar to zero point 
fluctuations in a Heisenberg magnet. The corresponding "Ising limit" where these fluc- 
tuations vanish is produced by adding an additional term to the Kondo interaction that 
couples the conduction electron isospin to the local moment as follows 



3 

The new term polarizes the composite order parameter, stabilizing a phase with A4 = 1. 
Decomposing the conduction electron into four Majorana components as in (15), then 




(16) 




(17) 



&j + Tj = —iiftj x Ipj, 



ffj - Tj = 2viptfj 



(18) 
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thus in the special case A = 1, the zeroth Majorana component explicitly decouples from 
the local moments, as in the mean field theory. 

One fascinating feature of the this isospin Kondo model, is its precise equivalence to 
the two channel Kondo model in the one impurity limit. In the one impurity version of 
this model the conduction sea is one dimensional and near the Fermi surface, spin-charge 
decoupling means that the isospin and spin of the conduction electron behave as two 
independent spin degrees of freedom, precisely emulating the two screening channels of the 
two-channel Kondo model. 18 This isomorphism is lost for multi-site or lattice models. In 
fact, NFL properties are far more stable in the isospin model, which manifestly preserves 
the decoupling of the neutral Majorana mode in the lattice at A = 1. Furthermore, the 
gap in the other three Majorana bands will preserve the decoupling of the neutral mode in 
a finite region about A = 1. Careful quantitative calculations of the zero point fluctuations 
are required to establish if the domain of attraction of the odd-u state extends all the 
way to the A = Kondo lattice (dimensionality plays an important role here). Clearly 
though, these simple considerations establish an important link between the existence of a 
lattice analog of the two channel Kondo fixed point and the development of odd-frequency 
pairing. 
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FIGURE CAPTIONS 



Fig. 1. (a) Scaling behavior of the two channel Kondo model; (b) illustration of a 
hypothetical new fixed point in the scaling trajectories of the Kondo lattice model. 

Fig. 2. C v /T and NMR 1/T\ calculated for a variety of chemical potential values illus- 
trating the failure to develop a Korringa relaxation in the presence of severe gaplessness. 
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